Pranck-Condon-Broadened Angle-Resolved Photoemission Spectra Predicted in 

LaMn0 3 
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The sudden photohole of least energy created in the photoemission process is a vibrationally excited 
state of a small polaron. Therefore the photoemission spectrum in LaMnOs is predicted to have 
multiple Franck-Condon vibrational sidebands. This generates an intrinsic line broadening » 0.5 eV. 
The photoemission spectral function has two peaks whose central energies disperse with band width 
« 1.2 eV. Signatures of these phenomena are predicted to appear in angle-resolved photoemission 
spectra. 
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The colossal magnetoresistance (CMR) effect [0 in 
doped manganese oxides has attracted a lot of atten- 
tion. The interplay of charge, orbital and magnetic or- 
der results in a very rich phase diagram S. The parent 
compound LaMnOs has orthorhombic symmetry at low 
temperature. The Mn +3 ion has d 4 {t2 g ,e g ) configura- 
tion with an "inert" t^g core (spin 3/2) and a half-filled 
doubly degenerate e s -type d orbital which is Jahn- Teller 
(JT) unstable. Ignoring rotation of the MnOg octahedra, 
which occurs below 1010 K, the JT symmetry breaking is 
cubic to tetragonal M at Tjt =750 K. The corresponding 
orbital order has x- and y-oriented E ff orbitals alter- 
nating in the x — y plane with wave vector Q = (tt, tt, 0). 
This in turn causes layered antiferromagnetic (AFA) or- 
der to set in at T N =140 K. 

The electronic structure of LaMnC>3 has been studied, 
for example, by photoemission [||J^] and by first princi- 
ples calculations 0. Still there is controversy about the 
nature of the low energy excitations, arising from the in- 
terplay between strong on-site Coulomb repulsion (which 
leads to magnetic order) and strong electron-phonon (e- 
p) interactions ||| (which lead to orbital order). 

When an electron is removed from the JT-ordered 
ground state, e-p coupling causes the hole to self-localize 
in an "anti-JT" small polaron state. In a previous pa- 
per H we have described the localized polaron in adi- 
abatic approximation. Residual non-adiabatic coupling 
allows the hole to disperse with band width narrowed 
by Huang-Rhys factor e ~ 3A / 4huJ » 10~ 4 . The photoe- 
mission process is sudden. The emitted electron with 
wavevector k leaves a hole in a lattice "frozen" in the 
unrelaxed JT state. Ignoring lattice relaxation, this hole 
would disperse with band width 2t sal eV (t is the hop- 
ping parameter), as shown on Fig. [l]. However, this is 
not a stationary state and must be regarded as a super- 
position of exponentially narrowed small polaron bands. 
Such bands have anti-JT oxygen distortions at each site, 
but a sufficient number of vibrational quanta are also 
excited such that the anti-JT distortion at time zero is 
"undone" . This is "Franck-Condon principle" . 



The measured spectrum at wavevector k will consist 
of a central 5-function at the energy of the frozen lattice 
(dispersive) band ei^(k), plus multiple vibrational side- 
bands at energy Si^(k) ± nhuj, with an overall Gaussian 
envelope whose width is approximately the polaron bind- 
ing energy. Franck-Condon broadening has been seen in 
photoemission spectra of solid nitrogen and oxygen |10|. 
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FIG. 1. Peak positions of the ARPES (a) at T=0 K and 
(b) at T=300 K along high- symmetry lines of the tetrag- 
onal Brillouin zone are shown by solid lines. Spin disorder 
at T above the Neel temperature Tn effectively reduces the 
hopping parameter t=(dder) by 2/3 and adds hopping in z di- 
rection which strongly affects dispersion. Error bars (FWHM 
of the photoemission peak) represent Franck-Condon broad- 
ening. Band structure (U=0 limit) peak positions Ai,2(fe) are 
shown by dashed lines. 

We find at each wavevector k the photoemission spec- 
trum has an intrinsic Franck-Condon broadening indi- 
cated by error bars in Fig. [l| The position of the maxi- 
mum disperses with /c-vector close to the "frozen" lattice 
spectrum. A qualitative picture of this process has been 
given by Sawatzky fill] in the context of high temper- 



1 



ature superconductors and by Dessau and Shen H for 
the manganites. The present paper gives an exact alge- 
braic prediction for the Angle-Resolved-Photoemission- 
Spectra (ARPES) of a model Hamiltonian for LaMnO,3. 

Our model Hamiltonian 0, first introduced by Millis 
Ji2| ], has hopping H e h electron-phonon Ti ep , and lattice 
Hl energies: 

Wei = tJ2 {[cl(i)c x {£± x)] + [x -> y] + [y^ z}} 

i,± 

Wcp = ~ 5 ^ ^,a(«i,a - Ut- a ) 

w L = J2( p U 2M + Ku U 2 )- (!) 

In these formulas c}JJt) creates a state with orbital 
ipa = |3ct 2 — r 2 >, where a = x,y,z. These three or- 
bitals span the two dimensional e g subspace and can 
be expressed in terms of the conventional orthogonal 
basis *2 = d x 2_ y 2, \& 3 = d 3z 2_ r 2 = ty z ] specifically 
= ±\/3/2*2 - *3/2. The resulting H c \ coincides 
with the nearest-neighbor two-center Slater-Koster fig] ] 
hopping Hamiltonian with overlap integral t — (dda) and 
(ddS) = 0. The hopping parameter t = 0.5 eV is chosen 
to agree with an ab initio e g band width of 1 eV . The 
e-p interaction 7i cp is modeled by a linear energy reduc- 
tion of an occupied i/j x orbital (fi£ x = c x (£)c x (£)) if the 
corresponding two oxygens in the ±x direction expand 
outwards, and similarly for y and z oxygens if i/i y or ip z 
orbitals are occupied. The strength of the e-p coupling g 
determines the JT splitting 2A = 1.9 eV, which is fitted 
to agree with the lowest optical conductivity peak |l4],|l5| . 
Static oxygen distortions 2u = ^2A/Afu 2 =0.296 A 
given by our model agree well with neutron diffraction 
data 0.271 A ||. For the lattice term Hl we use a sim- 
plified model where oxygen vibrations along Mn-O-Mn 
bonds are local Einstein oscillators. The displacement 
ui ta is measured from cubic perovskite position of the 
nearest oxygen in the d-direction to the Mn atom at £. 
The oxygen vibrational energy hui — Uy/ K/M — 0.075 
eV is taken from Raman data . In addition there is a 
large on-site Coulomb repulsion U and a large Hund en- 
ergy. These terms inhibit hopping except to empty sites 
where the t 2fl core spins are aligned correctly. 

In adiabatic approximation one can solve this problem 
for U=0 or U=oo. Both cases give a good description 
of the observed cooperative JT order. When U=0, the 
ground state wavefunction is: 



X\ 2 = A 2 + t 2 (2C 2 x + C x C y + 2C 2 ) ± t\C x + C y \ 



iGs ! o>=ni4 2 i vac 



> 
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A JT gap w 2A opens and the lower two bands of energy 
Ajr 1 , A^r 2 are filled. The photohole as initially created has 
energy: 



y jA 2 +4t 2 (C 2 -C x C v + C 2 ), (3) 

where C x _ y — cosfc^y and C z not entering at T=0 K. 
These bands are shown in Fig. |l](a) as dashed lines. 

At this point, at least in principle, could proceed nu- 
merically to find the polaronic energy lowering. However, 
it is both easier and more realistic to switch to U=oo. 
The ground state wavefunction: 



|GS,oo >= \\c x (£) JJ cj, (01 vac 



>, 



(4) 



has orbitals $x,y = —(^3 T ^ / 2)/v / 2 occupied singly on 
interpenetrating A and B sublattices. This is a fully cor- 
related state with zero double occupancy, while Eq. (||) is 
a band wavefunction in which two electrons are found on 
the same Mn atom with non-zero probability. For U> 6t 
« 3 eV, the state (||) has lower energy than state (|j) |j) . 
Neglecting creation of orbital defects with energy w 2A 
the "frozen" lattice approximation predicts photohole en- 
ergies A + £1,2 (k ), where £i )2 (fc) =±(t/2)(C x + C y ). We 
need to add a non-adiabatic treatment of the e-p cou- 
pling. The effective Hamiltonian Ti c s = TL C \ + Ti cp + Ti^ 
for the single hole is: 

Hi = J2{( d Ut±x)dx(£)+dUt±y)dx(£)) 

H% + Ht = X) 4 A + ^ Kct (a Q (£) + 

ai(£)-b a (£-a)-bt(£-a)) + ^>^K(4 (5) 

Here the operator d x (£) = cx(£) creates a hole in 
the JT ground state by destroying an electron on or- 
bital X at site £ (if £ £ A sublattice), and the oper- 
ator dy(£) — cy(£) creates a hole on B sublattice (if 
I G B). The phonon operators a^i) or W a (l) create 
vibrational quanta on the £ + a/2 oxygen, if £ G A 
or £ S B respectively. The e-p coupling constants are 
n XtV , x = X /A/12(1 + V3/2; 1 - v / 3/2; 1). The Hamilto- 
nian and all other energy parameters A and t in Eq. |E| 
are in units of hu>. The total Hamiltonian has an addi- 
tional term + H® p + which is obtained from Eq. 
@ by interchanging operators: 



d Y <-> dx,a x <-> by, a y «-> b x , a z «-> b z 



(6) 



and summing over the B sublattice. 

Following Cho and Toyozawa |H| we are able to diago- 
nalize Hamiltonian (|^) in a very large truncated basis of 
functions with a hole present on site £ and an arbitrary 
number of vibrational quanta p±z, P±y, P±z on the six 
displaced neighboring oxygens: 
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|^ A (£,M)>=4WII U ? Q (-^) 



Uh («a) 



(6t(^_ a ))P-c 



(4ffl) p+ ° 

|GS,oo > 



(7) 



The displacement operator U"(/c) = exp [— n(ae — a])] 
makes the H. cp + TL\^ part of the Hamiltonian diagonal. 
To get basis functions |4' B (£, {p}) > for holes on the 
B sublattice, the operators in Eq. ([?]) should be inter- 
changed according to Eq. (||). The next step is to build 
Bloch wavefunctions by Fourier transformation of the ba- 
sis functions Eq. (0). Then the Hamiltonian (||) will be 
diagonal with respect to fc-vector. The hopping term 
of the Hamiltonian H. c \ couples the \ 1 & A (£,{p}) > and 
|4' B (£',{p}) >-wavefunctions on the neighboring sites. 
The vibrational wavefunctions give a product of Huang- 
Rhys factors, but a shared oxygen contributes a non- 
factorizable overlap integral. However if one treats this 
shared oxygen as two independent atoms, one coupled to 
each site, then the Hamiltonian has a simple form: 



A 



-c+p'- 



P(p +a ,K a )P(p- a , K a )P(p' +a , K a )P(p'^ a , K a 



1/2 



. (8) 



where P(p,n) = exp(— K 2 )K 2p /pl is a Poisson distribu- 
tion. Since off-diagonal terms factorize, the analytical 
solution is available in this approximation: 



*i ,2 (£) 



{p}=0 a 

*£({p»±*f({p}) 



P(p +ai K a )P(P- a , K a ) 



1/2 



V2G'(£A)(Ea'(P+«' +P-af) - X\) 

The corresponding eigenvalues are: 

E{' 2 (fc) = j+ ^ 2 (fc), l + £l>2 (fc)G(4' 2 ) = 



(9) 



G{ Xx ) = e-^Y. 

p=0 



(3A/4)f 1 



pi 



p-x\ 



(10) 



The G'(x\) function in Eq. (||) is a first derivative of 
G(x\) and makes wavefunctions normalized. The correct 
solution needs a numerical diagonalization of the Hamil- 
tonian which explicitly includes vibrational states of the 
four shared oxygens. As can be seen on Fig. || the dif- 
ference is negligible between a typical spectrum obtained 
using approximation (|h]) and correct numerical treat- 
ment. The ground state of the Hamiltonian (H), with 
energy (from Eq. 10) E (k) = A/4 + a; (fc), corresponds 



to the anti-JT polaron. Its effective mass, deduced from 
d 2 xo(k) /dk 2 , provides a realistic alternative (exact for 
A — > or oo) to the available variational approaches E3] 
or exact quantum Monte Carlo simulations Jl9| . 
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FIG. 2. Energy dependence of the imaginary part of the 
Green function at T=0 and k = (7r/2, 0, 0). The solid line is 
numerical diagonalization of (^) in the subspace ( Q). The 
l ong -dashed line is the analytical approximation QIOJ) and 
(h_2|) . The spectral function consists of two asymmetric peaks 
with mean values A+ei,2. The zero-phonon line is seen at the 
adiabatic ground state polaron energy A/4 The dashed 
lines show the peak positions Ai^ of the uncorrelated electron 
theory Eq. (g). 

An ARPES experiment measures the spectral function 
A(k,oj) = —-^G(k,uj) with momentum k fully resolved, 
provided there is no dispersion in the direction perpen- 
dicular to the surface. Although LaMnOs is cubic, be- 
cause of the layered AFA magnetic structure, at low tem- 
peratures Mn e g electrons are two-dimensional and the 
spectrum can be measured: 



A(k,co) = J2\ </|d t (fc)|GS,oo> 
/ 



l 5(E — Ef). (11) 



The operator d'(k) excites a hole from the JT ground 
state. Summation over final eigenstates \f > includes 
summation over branch index i = 1,2 and number of 
phonons A = 0, l...oo. The sequence of delta functions in 
Eq. (|ll|) should be replaced by convolved local densities 
of phonon states, which we approximate by a Gaussian, 
5(E) -y exp(~E 2 /2-f 2 )/V2n"f- Substituting solution |j), 
(|l0|) into equation (O), we obtain the spectral function: 



A{k,uj) = 



G 2 {x\) 
G'(x x ) 



5{E-E\ 



(12) 



Equation @ along with @ gives the ARPES spec- 
trum normalized to J dwA(k, ui) = 1. The first energy 
moment of the spectrum ji7| coincides with the free hole 
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energy calculated in the "frozen" lattice approximation 
A + £1,2 (&) shown on Fig. |[ The edge of the spectrum 
corresponds to polaron creation at energy w A/4. This 
transition is weaker by 3 orders of magnitude than the 
peak at A + £i,2(fc)- 

At room temperature magnetic order is lost. The para- 
magnetic state is modeled by a mean field approxima- 
tion, namely scaling the effective hopping integral by 
2/3 and allowing hopping in ±z direction. This modi- 
fies the single particle energy band entering Eq. (^) to 
e 12 (fc) = t/3( - 2C Z ± (C x + C y )). But the JT orbital 
order is not destroyed at T=300 K and Franck-Condon 
broadening is still expected. When spins are disordered, 
k is not a good quantum number and additional broaden- 
ing is expected. Only phonon broadening of the ARPES 
along with peak positions are shown on Fig. |l|(b) . 

The angle-integrated spectrum, shown on Fig. I fol- 
low and high temperatures, has a width of about 1.2 eV 
and is almost temperature independent. The uncorre- 
cted (U=0) band structure, shown for comparison, is 
sensitive to magnetic order and therefore temperature 
dependent. 
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FIG. 3. Angle-integrated photoemission spectrum (solid 
line) at (a) T=0 K and (b) T=300K. At room temperature, 
spins are disordered, but JT order is not destroyed. For com- 
parison, the uncorrelated electron (U=0 limit) band structure 
density of states is shown by the dashed lines, with 2D Van 
Hove singularities at T=0 K. 

The existing photoemission data [^|j6) are consistent 
with our predictions. Higher resolution experiments are 
needed to test the theory and to unravel the nature of the 
lowest energy excitations in the LaMnC>3. To make such 
an experiment possible, a single domain sample (having 
2D dispersion at T=0K) is needed, with good control of 
oxygen concentration [E0| . 
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